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摘 要 : 基于 蒸 散发 (ET) . 归 一 化 植被 指数 (NDVI) 及 土地 利用 数据 利用 M-K 检 验 、 Sen 趋势 分 析 等 
方法 ,研究 2001 一 2015 年 黄河 流域 ET 时 空 分 布 及 不 同 植被 覆盖 /土地 利用 下 的 ET 变化 规律 。 结 果 
表明 :(1) 黄河 流域 年 均 ET 呈 东南 高 西北 低 的 空间 分 布 格局 ,与 植被 覆盖 和 土地 利用 的 关系 具有 
较 好 的 一 致 性 ;(2) 黄河 流域 ETNDVI 均 呈现 波动 式 增 长 状态 ,土地 利用 除草 地 外 其 他 类 型 面积 
均 有 增加 ,上 中 下 游 的 变化 具有 明显 差异 ;(3) 流域 大 部 分 地 区 ET 与 植被 呈正 相关 关系 ,植被 变化 
是 影响 ET 的 主要 因素 之 一 ;(4) 黄河 流域 土地 利用 类 型 对 ET 的 响应 时 空 差异 比较 明显 。 
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78 HU (Evapotranspiration ,ET) 主要 包括 地 表 
FEZ IS ZK TEA EZR , AT DAF AY Be wet 
域 需 水 和 植被 耗 水 情况 ” ,是 水 文 循 环 和 能 量 平衡 
的 关键 要 素 ”, 在 地 球 大 气 圈 , 水 轿 和 生物 疾 中 发 挥 
着 重大 作用 ”。 随 着 水 资源 短缺 和 生态 环境 问题 的 
研究 以 及 对 水 资源 合理 配置 的 迫切 需求 ,ET 的 时 空 
格局 变化 和 定量 分 析 显 得 至 关 重 要 “。 黄 河流 域 生 
态 环 境 脆 弱 、 水 土 流失 严重 、 水 资源 供需 矛盾 尖锐 
和 荒漠 化 问题 均 十 分 突出 ,因此 充分 认识 ET 和 
植被 覆盖 /土地 利用 的 时 空 变化 规律 及 其 响应 关系 ， 
不 仅 有 利于 理解 其 对 水 循环 的 影响 作用 机 制 ,也 可 
为 区 域 生 态 环境 问题 提供 科学 依据 ”。 

遥感 技术 是 开展 全 球 或 区 域 太 度 ET 与 植被 覆 
盖 / 土 地 利用 动态 监测 的 有 效 方法 " ,基于 遥感 的 
ET 研究 已 成 为 近年 来 该 领域 研究 热点 之 一 ”。 
国内 学 者 研究 了 博 斯 腾 湖 流域 ”高 原 地 区 “” RAE 
河源 区 ""、 黄 土 高 原 ”、 西 北 干旱 区 “等 区 域 ET 时 
空 分 布 规律 及 与 影响 因子 的 关系 ,并 取得 一 系列 的 
进展 和 突破 。Wang 等 从 中 发 现 低 植 被 覆盖 时 植 
被 的 变化 对 流域 ET 产生 负面 影响 , 随 着 植被 覆盖 度 
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的 增加 , 负 效 应 逐渐 向 正 效 应 转变 ;Wu 等 对 淮河 
流域 实际 ET 的 时 空 变化 特征 ,不同 覆盖 类 型 下 的 区 
域 实际 蒸 散 特 征 以 及 主要 影响 因子 进行 定量 分 析 ; 
Ma 等 2 采用 统计 学 方法 探讨 了 黄土 高 原 IT 时 空 变 
化 与 气候 因子 和 植被 变化 的 关系 ; 阮 宏 威 等 “研究 
了 1992 一 2015 年 中 亚 五 国土 地 和 覆盖 与 ET 时 空 变化 
特征 并 分 析 耕 地 蒸 散 耗 水 特征 ;Yan 等 3 研究 了 华 
北平 原 地 区 实际 ET 具有 明显 的 时 空 变 化 差异 且 植 
被 覆盖 度 .风速 和 气压 的 变化 是 影响 年 际 ET 趋 势 的 
主要 因素 ;Liu 等 乓 探讨 了 土地 利用 /覆盖 变化 和 和 气 
候 变 化 对 湿地 ET 的 影响 ,并 确定 ET 时 空 动 态 中 主 
影响 因子 的 重要 性 。 
本 文 基于 MODIS ET 归 一 化 植被 指数 (NDVI) 
和 土地 利用 数据 运用 M-K 检 验 .Sen 趋 势 分 析 等 方 
法 ,对 2001 一 2015 年 黄河 流域 及 上 中 下 游 ET 时空 
格局 变化 进行 分 析 和 评估 ,进一步 探讨 ET 与 不 同 植 
被 覆盖 /土地 类 型 之 间 的 关系 ,揭示 研究 区 内 ET 和 
植被 的 演变 规律 ,以 期 为 黄河 流域 生态 保护 和 高 质 
量 发 展 提供 科学 参考 。 
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河流 域 NDVI/A 


1 数据 与 方法 


1.1 研究 区 概况 

黄河 流域 位 于 我 国 干旱 、 半 干旱 和 半 湿 润 地 
区 , 介 于 96°~119°E、32°~42°N 之 间 , 流 域 面积 79.5x 
10' km"( 图 1)。 黄 河流 域 幅 员 辽 阔 , 地 势 西 高 东 低 ， 
河源 至 内 蒙古 河口 镇 为 上 游 地 区 , 河 段 落差 较 大 ; 
内 蒙古 河口 镇 至 河南 省 桃花 峪 为 中 游 地 区 ,地 形 为 
38 TER ,水 土 流失 严重 ;河南 省 桃花 峪 以 下 的 河 
段 为 下 游 地 区 , 主要 由 黄土 冲击 平原 形成 加 。 黄 河 
流域 属 大 陆 性 气候 ,从 东南 至 西北 依次 是 半 湿 润 气 


候 、 半 干旱 气候 和 干旱 气候。 全 年 气温 差异 明显 ， 
年 平均 降水 量 为 200~600 mm。 流 域内 地 貌 类 型 多 
ES) ,植被 呈 东 南 高 西北 低 的 空间 分 布 格局 , 主要 
土地 利用 类 型 为 林地 .草地 和 耕地 2? 。 


海拔 高 度 /m 
| 高 : 6128 


B io 


图 1 黄河 流域 位 置 示 意图 


Fig. 1 Location diagram of the Yellow River Basin 


1.2 数据 来 源 与 处 理 

本 文 使 用 的 MODIS 数据 来 源 于 美国 国家 航天 
Ja NASA 网 站 ,该 数据 集 特 征 参 数 完整 .时间 分 辩 率 
较 高 ,因而 在 区 域 时 空 分 布 特征 监测 方面 具有 一 定 
的 优势 。ET 数 据 来 源 于 MOD16A2 产 品 ,空间 分 辩 
率 为 500 m, 时 间 分 辩 率 为 8 d;NDVI 数 据 来 源 于 
MOD13Q1 产品 ,空间 分 辩 率 为 250 m, 时 间 分 辩 率 
为 16d; 土 地 利用 数据 来 源 于 MCD12Q1 产 品 ,IPGP 
全 球 植被 分 类 包括 17 种 土地 利用 类 型 。 本 文 分 析 
黄河 流域 主要 土地 利用 类 型 的 变化 特征 时 ,将 15 期 
土地 覆盖 数据 重 分 类 为 林地 、 草 地 、 永 久 湿地 ,农用 
地 建筑 用 地 和 水 体 ( 表 1)。 
1.3 研究 方法 
1.3.1 M-K 检 验 在 M-K 检 验 中 , 原 假设 HO 为 时 间 
序列 数据 (Xi, Xs, X) ,是 n 个 独立 的 、 随 机 变量 同 


ChinaXiv 合 作 期 刊 


上 地 利用 对 蒸 散发 时 空 变化 的 影响 


Rl 土地 利用 重 分 类 
Tab.1 Land-use reclassification 
重 分 类 IPGP 全 球 植被 分 类 
1 林地 1 常 绿 针 叶 林 
2 
3 落叶 针 叶 林 
4 落叶 阔 叶 林 
5 混交 林 
6 稠密 灌 从 
THEA, 
8 木 本 热带 稀 树 草原 
9 热带 稀 树 草原 
10 草 地 
16 fibi FEL UK 
3 湿地 11 永 久 湿地 
12 农 用 地 
14 农 用 地 /自然 植被 拼接 
5 建筑 用 地 13 城 市 和 建筑 区 
15 雪 和 冰 
17 水 


分 布 的 样本 ;假设 HI 是 双边 检验 ,对 于 所 有 的 k,j< 
n, 且 碾 j, 名 和 万 的 分 布 是 不 相同 的 ,检验 的 统计 量 5 
计算 如 下 式 : 


3 ; Sgn(X, - X,) (1) 
k=l j=k+l 
+1, (X,-X,)>0 
Sgn(X,-X,)=40, (X,-X,)=0 (2) 
-1, (X,-X,)<0 


5 为 正 态 分 布 ,其 均值 为 0, 方 差 Var(S)=nx(n- 
1)x(2n+5)/18。 当 n>10 时 ,标准 的 正 态 系统 变量 通 
过 下 式 计算 : 


Z -40, S=0 (3) 


v 
在 双边 的 趋势 检验 中 ,在 给 定 的 a 置信 水 平 上 ， 
如 果 1ZEZ wo, 即 在 置信 水 平 上 ,时 间 序 列 数据 存在 明 
显 的 上 升 或 下 降 趋 势 。 对 于 统计 量 Z,Z>0 时 是 上 升 趋 
势 ;Z<0 时 是 下 降 趋势 。|Zlz1.28、1.64 和 2.32 时 ,分 别 
表示 通过 了 信和 度 90% 9596 .99% 的 显著 性 检验 。 
13.2 Senda HA Me ”Sen 趋势 分 析 用 于 量化 变 
化 趋势 ,表示 单位 时 间 内 的 变化 量 ,计算 公式 为 : 
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$ 4 b 
s. (X-X 2[ ET 与 NDVI 土 地 利用 时 间 序 列 分 析 
B= median| 一 -一 (4) - 
i-j 2001—2015 年 黄河 流域 以 及 上 中 下 游 ET 和 


StH 1 <j<i<n, 当 B8>0 时 ,反映 了 时 间 序 列 呈 上 
升 趋势 ;反之 , 则 呈现 下 降 趋 势 。 


2 结果 与 分 析 


2.1 ET 与 NDVIL 土 地 利用 空间 特征 分 析 

黄河 流域 ET 与 NDVI 的 多 年 均值 空间 分 布 及 
2015 年 土地 利用 现状 如 图 2 所 示 。 高 荧 发 区 域 主要 
集中 在 黄河 源 区 .渭河 上 游 . 泾 河北 洛 河 、 延 河流 
域 等 ,而 风沙 区 的 燕 发 程度 较 低 。 黄 河流 域 不 同 区 
域 植被 空间 分 布 差异 明显 ,整体 上 呈现 由 东南 向 西 
北 逐 渐 减 小 的 变化 趋势 。NDVI 较 高 的 地 区 集中 在 
黄河 源 区 和 黄河 流域 中 下 游 , 风 沙 区 .河套 灌区 北 
部 ,十 大 孔 竞 等 地 区 植被 覆盖 较 低 。 黄 河流 域 主要 
的 土地 利用 类 型 是 草地 ,分 布 在 黄河 上 游 、 泾 河北 
部 清润 河北 部 以 及 汾 河西 北部 地 区 ; 泾 渭河 平原 
和 汾 河 平原 等 区 域 适合 农作物 的 生长 。 整 体 来 说 ， 
黄河 流域 ET 与 植被 覆盖 和 土地 利用 类 型 具有 一 定 
的 相关 性 。 


NDVI 年 际 变化 如 图 3 所 示 。 由 图 可 知 ,黄河 流域 
ET 整体 处 于 增长 阶段 ,2003 年 由 于 气温 高 于 常年 导 
致 ET 变化 剧烈 天 , 2004—2009 年 ET 增长 较为 平 
缓 ,2009 年 之 后 ET 增长 速度 明显 加 快 。 上 中 下 
游 ET 均 处 于 增长 水 平 ,黄河 流域 下 游 ET 的 波动 幅 
度 更 为 剧烈 。2001 一 2009 年 黄河 流域 下 游 蒸 发 量 
在 整个 研究 区 位 于 最 高 水 平 ,然后 依次 为 上 游 .中 
游 。 在 2009 年 之 后 黄河 流域 中 游 地 区 的 ET 增长 趋 
势 更 为 显著 ,在 2013 年 达到 377.35 mm ,黄河 上 游 地 
区 ET 处 于 研究 区 最 低 水 平 。 不 同 河 段 ET 虽然 差异 
较 大 ,但 年 际 波动 基本 一 致 。 黄 河流 域 的 植被 ND- 
VI 平均 值 常年 在 0.5 以 上 ,处 于 增长 水 平 。 年 NDVI 
最 小 值 . 最 大 值 出 现在 2001、2012 年 ,分 别 为 0.52、 
0.61。 黄 河流 域 下 游 地 区 植被 生长 旺盛 ,植被 指数 
均 在 0.75 以 上 ,变化 趋势 较为 平缓 ,无 明显 增长 趋 
势 ;2001 一 2015 年 黄河 中 游 地 区 植被 指数 处 于 平稳 
的 增长 状态 且 趋 势 明 显 ;黄河 流域 下 游 植被 指数 在 
0.45-0.55 之 间 ,NDVI 波 动 变化 情况 与 黄河 流域 整 
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Fig.2 Spatial distributions of ET, NDVI and land-use in the Yellow River Basin 
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图 3 黄河 流域 ET NDVI 时 间 序 列 变化 
Fig.3 Changes of ET and NDVI time series in the 


Yellow River Basin 


体 上 植被 变化 具有 较 好 的 一 致 性 。 总 体 来 说 ,黄河 
流域 不 同 河 段 植被 生长 变化 具有 明显 差异 ,ET 具有 
随 着 植被 盖 度 增加 而 增 大 的 趋势 思 。 黄 河流 域 ND- 
VI 和 ET 整体 上 呈 一 致 的 年 际 波动 ,各 个 河 段 的 ND- 
VI 和 ET 变化 趋势 大 致 基本 相同 ,在 一 定 程 度 上 表 
明年 ET 与 NDVI 具 有 较 好 的 相关 性 。 

2001 一 2015 年 黄河 流域 及 上 中 下 游 土地 利用 
类 型 面积 变化 情况 如 图 4、5 所 示 。 由 图 可 知 ,林地 、 
湿地 .农用 地 建筑 用 地 和 水 体面 积 持续 增 加 ,分别 
增加 了 0.8x10* km?, 0.03 x 10° km’, 1.81 X 10* km’, 
0.08X 10* km? .0.02x10: kn? ,而 仅 有 草地 面积 呈 下 降 
趋势 。 黄 河流 域 的 主要 土地 利用 类 型 具有 明显 的 
阶段 性 特征 ,不 同 土地 利用 类 型 之 间 的 变化 也 具有 
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图 4 黄河 流域 土地 利 
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明显 的 对 应 关系 。2001 一 2003 年 黄河 流域 农用 地 
面积 显著 增加 ,相应 时 期 的 草地 面积 显著 减少 ; 
2003 一 2011 年 黄河 流域 林地 、 农 用地、 建筑 用 地 面 
积 均 以 平稳 速率 上 升 ,同时 期 草地 的 面积 也 以 平稳 
的 速率 下 降 ;2012 一 2015 年 农用 地 水体 的 面积 处 
于 稳定 的 状态 ,该 时 间 段 内 草地 减少 的 趋势 逐渐 
WR. 

2001 一 2015 年 上 .中游 林 地 面积 的 增长 变化 速 
率 相对 稳定 ,而 下 游 地 区 的 林地 面积 波动 幅度 较 
大 , 且 在 2011 一 2013 年 林地 面积 减少 了 0.017x10' 
km’, Æ 2015 4 , ,林地 面积 基本 恢复 到 往年 最 高 水 
平 。 草 地 面积 一 直 处 于 减少 的 趋势 ,在 各 个 河 段 内 
变化 差异 较为 明显 。 上 游 地 区 的 草地 面积 由 2001 
年 的 44.91x10':km? 减 少 到 2015 4E AY 44.12x10* km’, 
面积 减少 速率 较 大 的 时 间 段 集中 在 2009 一 2013 年 ; 
中 游 地 区 草地 面积 在 2001—2012 年 均 以 平稳 速率 
减少 ,在 2012 年 以 后 基本 处 于 稳定 状态 ;2001 一 
2003 年 下 游 地 区 的 草地 面积 以 0.076x10' km?» a` AY 
速率 下 降 ,2003 年 以 后 草地 面积 呈现 波动 式 变化 ， 
无 明显 上 升 和 下 降 的 趋势 。 黄 河流 域 上 中 下 游 的 
湿地 面积 波动 情况 相对 较 小 , 均 以 缓慢 速率 增长 。 
黄河 流域 上 中 下 游 的 农用 地 、 建 筑 用 地 和 水 体 的 面 
积 处 于 稳定 的 变化 趋势。 
2.3 NDVI 变 化 对 蒸 散发 的 影响 分 析 

按 相 等 间隔 将 黄河 流域 年 均 NDVI 和 IT 分 为 $ 
级 来 分 析 NDVI 变 化 对 ET 的 影响 ( 表 2、 表 3)。 研 究 
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Fig.4 Area changes of land-use types in the Yellow River Basin 
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Fig. 5 Area changes of land-use types in upper, middle and lower reaches of the Yellow River Basin 
表 2 研究 区 各 河 段 及 全 区 平均 NDVI 及 NDVI 等 级 比例 
Tab.2 Average NDVI and NDVI grade ratio of each river section and the whole area in the study area 
. m NDVI 等 级 /% 
河 段 平均 NDVI a E 
<0.2 0.2~0.4 0.4~0.6 0.6~0.8 >0.8 
上 游 0.53 7.04 25.57 25.54 30.61 11.16 
中 游 0.69 0.13 4.71 21.59 47.47 26.10 
下 游 0.76 0.11 1.07 7.96 48.47 42.33 
全 区 0.59 4.62 18.02 23.31 36.57 17.14 
HE: NDVI, 归 一 化 植被 指数 。 
表 3 研究 区 各 河 段 及 全 区 平均 ET 及 ET 等 级 比例 
Tab.3 Average ET and ET grade ratio of each river section and the whole area in the study area 
ET 等 级 /4 
河 段 平均 ET/mm Eau 
«200 mm 200-400 mm 400-600 mm 600-800 mm 2800 mm 
上 游 318.70 30.41 32.99 36.51 0.09 0 
中 游 316.69 13.69 63.19 22.45 0.71 0 
下 游 337.70 4.60 74.70 20.14 0.56 0 
全 区 319.24 23.99 44.30 31.42 0.30 0 


区 多 年 平均 NDVI 为 0.59,22.64% 的 区 域 NDVI<0.4， 
仅 有 17.14% 的 区 域 NDVI>0.8,NDVI 介 于 0.4~0.8 的 
区 域 占 总 流域 的 59.88%。 上 中、 下 游 NDVI>0.6 的 
区 域 所 占 比例 越 来 越 高 ,表明 从 东 至 西 植被 覆盖 度 


VE: ET, ALE 


越 来 越 高 。 下 游 地 区 的 NDVI 达 到 最 高 为 0.76, 该 
河 段 超过 90% 的 区 域 NDVI>0.6。 黄 河流 域 年 均 ET 
为 319.24 mm,44.30% 的 区 域 ET 在 200~400 mm 之 
间 ,ET 在 0~200 mm 和 400~600 mm 分 别 占 流域 的 
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23.99% , 31.4290 , 而 仅 有 0.30% 的 区 域 ET 值 达到 
600 mm 以 上 。ET 处 于 200~400 mm 阶段 的 面积 所 
占 比例 随 着 上 游 中游、 下 游 逐 渐 增 加 ,表明 ET 程度 
从 东 至 西 越 来 越剧 烈 。 下 游 地 区 95.40% 的 区 域 ET 
大 于 200 mm ,年 均 ET 最 大 。 黄 河流 域 下 游 是 NDVI 
和 ET 高 值 的 集中 分 布 区 ,中 低 程 度 植被 覆盖 和 ET 
分 布 在 流域 中 上 游 。 总 体 来 说 ,研究 区 NDVI 与 ET 
的 空间 分 布 基本 一 致 ,同时 表明 了 ET 的 大 小 与 植被 
覆盖 的 高 低 具有 较 好 的 相关 性 。 

根据 M-K 检 验 与 Sen 趋 势 分 析 的 方法 对 研究 区 
NDVI 和 ET 逐 像 元 计算 , RAG NDVI 5j ET HY M-K $È 
验 统计 量 Z. 值 及 年 际 变 化 速率 6 值 结 果 , 并 对 NDVI 
和 ET 的 变化 趋势 及 变化 速率 进行 分 类 和 划 等 ,统计 
其 比例 特征 ( 表 4)。 人 研究 发 现 : 中 游 的 NDVI 和 ET 
的 变化 速率 最 大 ,分 别 为 0.007.8.26, 呈 显著 增加 的 
变化 趋势 ,通过 a=0.05 显著 检验 ;上 游 地 区 植被 指 
数 的 之 值 为 1.59, 达 到 显著 化 增加 水 平 ,而 ET 变化 
是 不 显著 水 平 增加 ;下 游 地 区 NDVI 的 变化 速率 仅 
为 -0.001, 生 长 情况 处 于 稳定 状态 ,ET 变化 则 处 于 
平稳 增加 阶段 。 从 流域 整体 来 看 ,植被 NDVI 和 ET 
变化 情况 呈现 不 显著 增加 趋势 。 

由 表 5、 表 6 可知, 黄河 流域 NDVI 变 化 率 主 要 
介 于 0~0.01 之 间 , 其 比例 为 69.54% ,变化 率 小 于 0 
和 变化 率 大 于 0.01 的 比例 分 别 为 15.32% 和 
15.14%。 上 游 地 区 NDVI 变 化 速率 等 级 在 -0.005~0 
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之 间 占 其 区 域 的 51.55%, 人 处 于 缓慢 增长 阶段 ;中 游 
地 区 植被 增长 趋势 明显 ,30.41% 的 区 域 NDVI 增 长 
速率 大 于 0.01; 下 游 47.96% 区 域 的 植被 呈 缓 慢 减少 
的 趋势 。 黄 河流 域 ET 变化 率 主要 介 于 0~12 mma", 
其 比例 为 72.57% ,变化 率 小 于 0 mm:a ' 及 变化 率 大 
于 12 mm*a! 的 比例 分 别 为 14.93% 和 12.49%。 中 游 
是 NDVI 及 ET 增加 速率 最 快 的 集中 区 域 ,其 NDVI 变 
化 速率 大 于 0.005 和 ET 变化 速率 大 于 4mma' 比 例 
分 别 达 到 了 63.02% 和 77.95%。 上 游 和 下 游 地 区 ND- 
VI 与 ET 的 主要 呈现 增加 趋势 ,但 变化 相对 幅度 较 小 。 

利用 2001 一 2015 年 黄河 流域 ET 和 NDVI 的 长 
时 间 序 列 数据 分 析 其 相关 系数 (图 6)。 由 图 6 可 知 ， 
黄河 流域 大 部 分 地 区 ET 与 植被 均 呈 正 相关 关 系 , 相 
关 性 较 好 的 地 区 主要 集中 在 河 龙 区 间 黄 丘 区 、 泾 渭 
河上 游 .内 蒙古 黄河 段 东北 部 等 区 域 。 黄 河流 域 中 
游 南 部 地 区 降雨 量 较 多 ,植被 生长 茂盛 ,蒸发 程度 
较 高 ,ET 与 NDVI 变 化 的 趋势 较为 一 致 。 相 关系 数 
相对 较 小 的 区 域 主要 集中 在 黄河 源 区 ,该 区 域 海 拔 
高 地 貌 多 系 山岭 及 草地 高 原 ,植被 覆盖 度 较 低 , 但 
是 ET 程度 相对 较 高 。 农 用 地 区 域 的 ET 和 NDVI 相 
关 性 较 低 。 整 体 来 说 ,ET 的 变化 情况 与 植被 生长 状 
态 具 有 一 定 的 相关 性 ,植被 的 长 势 是 影响 ET 的 主要 
因素 之 一 。 
2.4 土地 利用 变化 对 蒸 散 发 的 影响 分 析 

图 7 为 黄河 流域 各 河 段 土地 利用 和 ET 的 变化 


表 4 研究 区 全 区 及 不 同 河 段 NDVI 和 ET 变化 趋势 


Tab.4 Variation trends of NDVI and ET in the whole study area and different river sections 


IE NDVI ET 
Z. B 趋势 Z. B 趋势 

上 游 1.59' 0.004 显著 增加 0.88 4.02 不 显著 增加 

中 游 2.10 0.007 显著 增加 1.65; 8.26 显著 增加 

下 游 0.06 -0.001 不 显著 减少 0.68 3.55 不 显著 增加 

全 区 0.58 0.005 不 显著 增加 1.09 1.87 不 显著 增加 


注 :Z、p 分 别 表 示 M-K 检 验 统计 量 、 年 际 变化 速率 ;* 代 表达 到 显著 水 平 , P<0.05。 


表 5 研究 区 全 区 及 各 河 段 NDVI 不 同 变 化 速率 等 级 比例 


Tab.5 Proportions of NDVI with different rate of change in the whole area and river sections of the study area 


. NDVI 变 化 速率 等 级 /% 
河 段 

<-0.005 -0.005~0 0~0.005 0.005~0.01 >0.01 
上 游 0.81 13.27 51.55 24.48 9.89 
中 游 2.98 6.84 27.15 32.61 30.41 
下 游 14.77 33.19 36.78 11.96 0.3 


全 区 2.3 13.02 43.62 25.92 15.14 
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表 6 研究 区 全 区 及 各 河 段 ET 不 同 变化 速率 等 级 比例 
Tab. 6 Proportions of ET with different rate of change in the whole area and river sections of the study area 
ET 变化 速率 等 级 /mm ‘a! 


me -4 -4-0 0-4 4-8 8-12 >12 

上 游 1.79 16.44 39.05 23.82 11.90 7.00 
中 游 2.53 3.02 16.05 26.83 25.50 26.07 
下 游 9.27 15.50 24.01 28.27 16.41 6.53 
全 区 2.46 12.47 31.51 24.97 16.09 12.49 


特征 。 从 流域 整体 来 看 ,草地 作为 研究 区 主要 的 土 
地 利用 类 型 随时 间 变 化 面积 逐渐 减少 ,其 他 土地 利 
用 类 型 面积 则 逐渐 增加 ,但 6 种 土地 利用 类 型 的 蒜 
发 量 均 呈 增加 的 趋势 。 黄 河流 域 70% 以 上 的 草地 
位 于 研究 区 的 上 游 ,各 类 土地 利用 类 型 面积 的 变化 
与 整个 流域 较为 一 致 ,同样 蒸 发 量 也 随 着 时 间 变 化 E ciis 
逐渐 增 大 ,研究 区 内 的 林地 和 农用 地 主要 分 布 在 黄 mm 


To] fce DC 6 BL FA FES RI HIZR ERI Oana ane ene 
间 变 化 基本 一 致 。 黄 河流 域 下 游 土 地 利用 类 型 以 Fig.6 Spatial distribution of correlation coefficients 


between ET and NDVI 


农用 地 为 主 ,各 类 土地 类 型 的 蒸发 程度 均 大 于 中 上 
游 , 且 湿 地 和 林地 的 蒸发 量 较 大 。 典 型 流域 中 浑 ” 河 下 游 主 要 的 土地 利用 类 型 是 农用 地 ,清润 河 、 延 
河 \ 十 大 和 孔 竞 .皇甫 川 窜 野 河 、 秃 尾 河 、 无 定 河 均 以 ” 河流 域 以 林地 为 主 , 藻 发 量 从 西北 向 东南 逐渐 增 
草地 为 主要 的 土地 利用 类 型 , 泾 河北 洛 河 以 及 渭 。 加 。 综 上 所 述 , 土 地 利用 是 影响 ET 的 重要 因素 ”， 
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Fig. 7 Land-use types and evapotranspiration changes in the Yellow River Basin 
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黄河 流域 土地 利用 类 型 对 ET 的 响应 时 空 差异 比较 
明显 。 明 确 ET 实 际 时 空格 局 变化 特征 ,对 于 合理 分 
配水 资源 .提高 水 资源 利用 率 具有 重要 意义 。 


3 结论 


目前 ,对 于 ET 的 研究 较 多 利用 MODIS 产品 或 
者 遥感 估算 模型 来 分 析 ET 的 时 空 变 化 特征 , 较 多 集 
中 在 小 尺度 范围 ,而 对 黄河 流域 iT 的 研究 较 少 。 
MODIS ET 数据 集 在 充分 考虑 植被 蒸腾 oe E ER A 
蒸发 和 土壤 水 分 蒸发 的 基础 上 能 够 满足 区 域 尺度 
的 研究 。 本 文采 用 MOD16A2 产品 可 以 有 效 地 反映 
黄河 流域 ET 的 变化 状况 ,有 助 于 促进 区 域 生态 环境 
的 保护 和 水 资源 合理 配置 。 基 于 MODIS 的 ET、ND- 
VI 和 土地 利用 数据 研究 了 2001 一 2015 年 黄河 流域 
ET 时 空 变化 特征 ,采用 趋势 分 析 和 相关 分 析 探 索 了 
上 中 下 游 ET 与 NDVI 响 应 关系 ,统计 分 析 了 黄河 流 
域 不 同 土地 利用 类 型 的 ET 变化 规律 。 主 要 结论 
如 下 : 

(1) 2001—2015 年 黄河 流域 及 上 中 下 游 ET 均 
处 于 上 升水 平 ,下 游 ET 的 波动 幅度 更 为 剧烈 。 人 研究 
区 年 均 NDVI 大 于 0.5 且 植被 变化 呈 波 动 式 增长 状 
态 , 不 同 河 段 植被 年 生长 变化 差异 较为 明显 。 

(2) 黄河 流域 NDVI 与 ET 的 空间 分 布 基本 一 
致 , 在 中 游 地 区 ET、NDVI 均 处 于 显著 增加 的 水 平 。 
大 部 分 地 区 ET 与 植被 均 明 正 相关 关系 ,植被 的 变化 
是 影响 ET 的 主要 因素 之 一 。 

(3) 研究 区 在 2001—2015 年 除草 地 以 外 其 余 
土地 利用 类 型 面积 均 呈 现 增加 的 趋势 ,量具 有 明显 
的 阶段 性 特征 ,不 同 河 段 的 土地 利用 类 型 面积 的 变 
化 趋势 均 有 一 定 的 差异 ,整体 来 说 黄河 流域 土地 利 
用 向 良好 的 方向 发 展 。 

(4) 不 同 河 段 土地 利用 类 型 对 于 ET 的 响应 具 
有 一 定 的 差异 ,黄河 流域 下 游 地 区 莹 发 量 明 显 
较 大 。 
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PERSE: 黄河 流域 NDVI 土 地 利用 对 蒜 散 发 时 空 变 化 的 影响 


Effects of NDVI/land-use on spatiotemporal changes of 
evapotranspiration in the Yellow River Basin 


GU Jiahe", XUE Huazhu, DONG Guotao”, ZHOU Lijuan’, 
LI Jingru', DANG Suzhem, LI Shangzhi' 


(1. School of Surveying and Land Information Engineering, Henan Polytechnic University, Jiaozuo 454000, Henan, China; 
2. Yellow River Institute of Hydraulic Research, Zhengzhou 450003, Henan, China; 3. Heihe Water Resources and Ecological 


Protection Research Center, Lanzhou 730030, Gansu, China) 


Abstract: The Yellow River Basin is mainly distributed in arid, semi-arid, and semi- humid areas with a fragile 
ecological environment, which has become the most serious soil erosion area in China and even in the world. 
Evapotranspiration (ET) is an essential part of the land water and energy cycle, playing a vital role in the global 
ecosystem regulation and the hydrological cycle. ET is mainly composed of vegetation and transpiration of soil 
water and different vegetation cover or land-use types, showing different spatial-temporal distribution characteristics. 
Accurate exploration of the ET’ s spatial-temporal distribution and its response relationship is conducive to mastering 
regional surface water and heat balance law. Remote sensing has been widely used in ET’ s dynamic monitoring 
on a global or regional scale. Analyzing the ET’ s spatial-temporal changes in the Yellow River Basin promotes 
understanding the impact of vegetation and land-use changes on the water cycle and the rational allocation of water 
resources. On the basis of ET- MODIS data, normalized difference vegetation index (NDVI), and land- use 
products as data sources, we studied the ET’ s spatial-temporal change pattern in the Yellow River Basin from 
2001 to 2015 using the Manner- Kendall test and Sen' s trend analysis. We discussed ET' s impact and the 
different change characteristics under different NDVI conditions and land-use types. We found the following 
results: (1) The ET distribution pattern in the Yellow River Basin is obvious, gradually decreasing from southeast 
to northwest. The average ET value in the study area is 319.24 mm, and the average NDVI is 0.59. The high 
values of ET and NDVI are concentrated on the downstream areas, and their spatial distribution is in good 
consistent. (2) The time series of ET and NDVI in the Yellow River Basin all showed a fluctuant growth state, 
and the changes in each river segment were significantly different from 2001 to 2015. In inter-annual changes, the 
trends of ET and NDVI are the same, and they have a significantly increased level in the middle reaches. (3) 
Vegetation change is the main factor affecting ET, positively correlating with ET in most Yellow River Basin 
areas. (4) From 2001 to 2015, the area of land cover types in the Yellow River Basin except the weeding areas 
showed an increasing trend, in which the ET response in different river reaches was different to some extent. The 
spatial-temporal differences in land-use types in response to ET are obvious. It is of great significance for ecosystem 
protection and water resources management to comprehensively analyze the variation law of evaporation in the 
upper, middle, and lower reaches of the Yellow River Basin and provide a reference for promoting the sustainable 
development of the regional economy and society. 


Key words: evapotranspiration; vegetation cover; NDVI; land use; Yellow River Basin 


